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Recently, based on the refined crystal structure of Pro.6Cao.4MnC>3 from neutron diffraction, 
' Daoud-Aladine et al.[Phys. Rev. Lett. 89, 97205(2002)] have proposed a new ground state structure 

for the half-doped manganites Ro.sCao.sMnOa, where R is a trivalent ion like Bi, La, Pr, Sm or Y. 
Their proposal describes the CE magnetic structure attributed to these materials as an arrangement 
s ! ' of dimers along the ferromagnetic Mn zig-zag chains that form it. The formation of dimers is in 

!_i ■ accordance with previous theoretical calculations based on simplified models for these compounds. 

However, the dimers' proposal is in conflict with the Goodenough-Kanamori- Anderson rules, which 
give a coherent description of many transition metal insulating compounds and predict the coex- 
istence of Mn 3+ and Mn 4+ ions in equal parts in the half-doped manganites. On the other hand, 
Rivadulla et al. [Phys.Rev.B 66,174432(2002)] have studied several single crystal samples of half- 
doped manganites and propose a phase diagram in terms of the tolerance factor which contains 
both types of structures. In the present work we have calculated the magnon dispersion relations 
for the CE magnetic structure, arising for each type of proposal: the charge ordered and the dimer 
phases, respectively. We consider a three-dimensional unit cell containing 16 spins, and compare the 
• magnetic excitations along different paths in the first Brillouin zone. We conclude that measure- 

ment of the magnon dispersion relations should allow a clear distinction between the two proposals, 
■4_> , predicting qualitative differences arising along specific directions of propagation in the first Brillouin 

zone. 
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Half-doped manganites of stoichiometry Ri-^Ca^MnOa, x= 1/2, where R is a trivalent ion like Bi, La, Pr, Sm 
or Y, display a rich phase diagram depending on their composition 1 : in the absence of magnetic field, charge- 
ordered antiferromagnetic insulator phases (COA), ferromagnetic metallic (FM), antiferromagnetic insulator, and 
paramagnetic phases have been reported. Strong colossal magnetoresistance appears by driving the materials from 
COA to FM by application of moderate external magnetic fields. 

Below the Neel temperature the half-doped compounds are found in the CE magnetic structure 2 , in which ferromag- 
netic Mn zig-zag chains contained in parallel planes are found in an antiferromagnetic arrangement. The description of 
the CE magnetic structure found in most half-doped manganites is based on the well known Goodenough-Kanamori- 
Anderson (GKA) rules 3 , implying a checkerboard ordering of Mn 3+ /Mn 4+ ions in the crystal 4 , as shown in Fig.l. 
■ ( +orbital) ordered (CO) scenario for half-doped manganites in 
£f) ' However, Daoud-Aladine et al. 5 have shown recently that the crystal structures of Pro.6oCao.4oMn03 single crystals 
obtained from neutron diffraction data can not be interpreted in terms of Mn 3+ /Mn 4+ charge ordering. Rather, the 
data suggest an intermediate valence of 3.5 for all Mn ions. Furthermore, the XANES results on the Mn K-edge of 
Lao.soCao.soMnOs by J.Garcia et al. 6 suggest also the same valence for Mn. On the basis of these results, Daoud- 
Aladine et al. propose that the zig-zag ferromagnetic chains contained in the CE magnetic structure are formed by 
a succession of Mn dimers (or Zener polarons) 5 , such as indicated in Fig. 2. Mn ions forming the dimers share one 
electron and point in the same direction forming a large spin unit. The chains order antiferromagnetically between 
them within the Mn structure. Theoretically, such situation had been studied previously in one 7 , two 8 and three 
dimensional 9 structures as a consequence of the competition between double exchange and superexchange. 

More recently, Rivadulla et al. 10 have proposed a new phase diagram for the half-doped manganites in terms of 
• i-H . the tolerance factor t. In the new phase diagram both phases appear: the charge ordered phase is found for t below 
1 a critical value t c = 0.975, and the dimer phase for t larger than t c . Both (insulating) phases are separated by a 
unique (metallic) ferromagnetic phase. Contrary to the usual metallic ferromagnetic phase of hole doped manganites 
which is stabilized by pressure, the ferromagnetic phase that separates the two insulating phases at half filling can be 
suppresed by pressure. 

Since both insulating phases proposed for the CE magnetic structure, the charge ordered and the dimer one, 
have some common characteristics like: their insulating character, both are metamagnetic going into a metallic 
ferromagnetic phase under moderate magnetic fields, include zig-zag ferromagnetic chains, etc., we decided to calculate 
their magnetic excitation spectra as a mean to distinguish between them, and also as a way to check the difference of 
the magnetic interaction along the chains in the proposed models. 

Therefore, we have considered a basic three-dimensional unit cell containing 16 spins (Mn), half of them in each of 
two consecutive planes, as shown in Fig. 3. For the charge ordered phase, two different values of spin (representing 



1 



the Si = 2 and S 2 = 3/2 spins of Mn 3+ and Mn 4+ , respectively) would appear, in checkerboard arrangement like 
shown in Fig. 1. The interactions are superexchange couplings between all Mn ions: one ferromagnetic coupling along 
nearest neighbour spins on a chain (F'=F), and two antiferromagnetic nearest neighbour couplings: one inter-chain 
coupling (A) and one inter-plane (A'). On the other hand, for the dimer phase, only one value of spin (S\ = S 2 ) would 
be present, while the ferromagnetic coupling between nearest neighbour spins along a zig-zag chain now would take 
two different values: the interaction between the spins forming a dimer is double exchange (F'), while a superexchange 
interaction (F) couples nearest neighbour spins in different dimcrs. 

This state of things can be described by a model composed of a series of appropriate nearest neighbour Heisenberg 
coupling terms. With the notation for magnetic couplings introduced in Fig. 3, the following Hamiltonian has been 
considered: 
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We have identified the position of the spins by indices: a = k) to locate the 3D unit cell in the Bravais lattice, 
and site labels n(= 1, .., 16) to distinguish the spins inside one unit cell. Capital letters D denote a dimer, NN 
nearest neighbour sites, C one chain, and P one plane. F, F' , A and A' are positive. For the charge-ordered phase 
the ferromagnetic couplings along a zig-zag chain are F'=F, while each coupling term will include spins of different 
magnitude (due to the checkerboard arrangement shown in Fig. 1). 

To calculate the magnetic excitations at low temperatures we have used the Holstein-Primakoff transformation 
in linear spin wave approximation 11 , which in our case demanded the introduction of 16 different kinds of boson 
operators, one for each spin inside the unit cell of Fig. 3, distinguishing the spin-up and spin-down sublattices. That 
is, we introduced for the spin-up sublatticc: 
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where n = 1,2,3,4,13,14,15,16 with the notation of Fig. 3, and S n is the magnitude of spin n. Without loss of 
generality, we can put: S n = Si, for n = 1, 3, 14, 16, and S n = S2, for n = 2, 4, 13, 15 ( Si = S2, for the dimer phase). 
Similarly for the spin-down sublattice, we introduced: 
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C- — 



/2S n b^ na 
-S n + bi_ n b, 



n,a u n,a , 
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where n = 5,6,7,8,9,10,11,12 with the notation of Fig. 3. Without loss of generality, we put: S n = Si, for n 
5,7,10, 12, and S n = S 2 , for n = 6,8,9,11. 



Introducing the Fourier transform of the boson operators (which for compactness we will denote here: a - = 
CLif 1 ; = bn^) , we then determined the Hamiltonian, which for the dimer phase (Si — t. 

n,k 

form: 



71, fc 

S) takes the following 



H=[F' + F + 2(A + A')} [-8S*N + S (E^i-4,13-16 4«n + E fc - n=5 _ 12 
- SF 'Ek [4«3 + b\b 6 + b f m bn + a\ 3 a u + e lk * a *a\a 4 + e ife ^al 5 a 16 + e - lk y a yb% + e~ lk * a *b\b 12 + H.c 
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+ SAJ^k [aib 6 + a 3 b 6 + a 3 b 8 + b 7 a 4 + b 9 a 14 + bna 14 + bna 16 + ai 5 b 12 + e~ lk * a *aib 8 + e~ lk * a * a 13 b 12 + 
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where N denotes the total number of 16-spin unit cells as in Fig. 3. 
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A similar expression can be written for the Hamiltonian of the charge-ordered phase, taking into account the two 
values of spin which will appear (Si ^ S2) and F = F' for that case. 

The magnon excitations of the system for the dimcr and the CO phases can then be obtained by diagonalization 
of the corresponding Hamiltonian. Concretely, we can obtain them by diagonalization of the following 16x16 general 
matrix (multiplied by S = V^iS^) , which for appropriate values of its parameters (as mentioned above) describes the 
two phases as particular cases: 
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where the matrix columns (or rows) are related to the 16 boson operators, arranged in increasing order of their 
site label, and we have used the following notation for the matrix elements: Ai = S2[F + F' + 2(A + A')]/S; 
A 2 =S 1 [F + F' + 2(A + A')]/S ; $ Q - cxp (ik a a a ) , (a = x,y) ; 9 X = (1 + exp (ik z a z )). 

In Fig. 4 we exhibit our result for the 16 magnon excitation branches of the charge ordered (CO) phase, plotted along 
8 different paths throughout the first Brillouin zone, using a typical set of relative values for the exchange coupling 
parameters for the system (in the present work we take the interdimer exchange coupling F as the energy unit, setting 
the overall energy scale) and a ratio of spin values S1/S2 = 1-3, such as that for Mn 3+ /Mn 4+ magnetic moments. 
The trajectories along the first Brillouin zone depicted extend between the following points: T = = (0,0,0), 
X = (7r/a x ,0,0), Y = (0, ir/a y , 0), Z = (0, 0, n/a z ), M = (-n/a x , n/a y , 0) and L = (ir/a x , Tr/a y , n/a z ). 

In Fig. 5 we exhibit the results calculated for the magnons corresponding to the dimer phase, along the same first 
Brillouin zone trajectories. We have now considered: S1/S2 — 1, F' — 2F, and the rest of the parameters unchanged 
for the sake of simplifying the comparison with the CO magnons. We have verified that the same magnon excitations 
are obtained for this phase, if one chooses a different relative array of the dimers on neighbour zig-zag chains. 

Comparison of Figs. 4 and 5 reveals that a qualitative difference in the magnon branches predicted for the CO 
and the dimer phases can be clearly pointed out. Namely, that along the T — ► Y and especially along the Y — ► M 
paths in the first Brillouin zone 8 CO magnon branches are doubly degenerate, while additional degeneracies are 
predicted for the respective magnons in the dimcr phase. For A' < A six distinct magnon branches (4 of them doubly 
degenerate, but other 2 with degeneration four) appear along T — > Y (while for A = A' as in Fig. 5 there is a slight 
lifting of the extra degeneracy, being 8 branches observable) , while 4 magnon branches (with degeneration four) 
appear along Y — ► M for the dimer phase. This splitting or lifting of the dimcr phase degeneracies observable in 
the CO phase increases proportionally with the magnitude of the antiferromagnctic in-planc A coupling (the effect 
disappears for A = 0), being almost unaffected by the values of the other magnetic couplings. 

The anisotropy exhibited by the magnon dispersion relations (of Figs. 4 and 5) along T — ► X with respect to those 
along r — > Y can be understood considering the spin ordering along the x and y directions (see Fig. 3). Along x, 
the direction of the zig-zag chains, the system can be viewed as composed by consecutive chains: one of parallel spins 
(e.g. up), followed by another with alternating spins, then again one of parallel spins (down), etc. While along y the 
system has more symmetry: it can be viewed as consecutive chains, all of them formed by alternating spins. 

Apart from the qualitative differences between the magnon spectra of the CO and dimer phases which were pointed 
out above (mostly related to the magnitude of the antiferromagnetic in-planc coupling A), there are also quantitative 
variations of the spectra related to the values of the other magnetic couplings considered. In particular, in Fig. 6 we 
exhibit the magnon dispersion curves for a set of parameters more suitable for the description of the dimer phase: with 
an intradimer ferromagnetic coupling (F 1 ) an order of magnitude larger than the interdimer ferromagnetic coupling 
(F). This corresponds to the large estimated values (t <~ 100 mcV) for the effective hopping between Mn ions in 
these compounds 12,13 . Recent estimations of the ferromagnetic interdimer exchange coupling F in these compounds 
consider it to be of the order of 10 mcV 13 . For strong intradimer ferromagnetic coupling as in Fig. 6, a large energy gap 
would then separate the magnons between higher (lying approximately around O.leV) and lower (around 0-0.01eV) 
energy modes, being the latter possibly the only ones easier to probe experimentally. 

Summarizing, we have calculated the magnetic excitations to be expected for the charge ordered and the dimcr 
phases which have been proposed for the CE magnetic structure of the half-doped manganites. We predict qualitative 
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differences between the magnons related to specific exchange couplings, which should allow to distinguish between 
these two phases in neutron scattering experiments. Not only are there lifted degeneracies in the CO phase along 
certain first Brillouin zone paths with respect to the dimer phase, but the presence of dimers with strong intradimcr 
exchange coupling should result in a smaller number of magnons being detected by the experiments at low energies, 
while a larger number of magnetic excitations should be easily detected if the charge-ordered phase was the one being 
measured. 
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Fig. 1. Schematic picture of the Mn 3+ /Mn 4+ charge ( +orbital) ordered (CO) scenario for half-doped manganites in the CE 
magnetic structure. Large arrows denote S = 2 ( Mn 3+ ) spins, while small arrows S = 1.5 ( Mn 4+ ) spins. I and II are coplanar 
ferromagnetic zig-zag chains, chain III is directly below I on a consecutive plane. 

Fig. 2. Schematic picture of the dimer phase for half-doped manganites in the CE structure. All spins correspond to Mn 
ions with an intermediate valence of ~3.5. Thick dashed lines indicate dimers. 

Fig. 3. Schematic plot of the three-dimensional 16-spin unit cell used ( Bravais lattice parameters: a x — a y — v^ao, a z = 2c). 
Magnetic couplings between nearest neighbour spins: ferromagnetic "intradimer" F' (thick dashes) and "interdimer" F (solid 
line) along the zig-zag chains (F'=F, in CO phase), antiferromagnetic interchain A for coplanar spins (e.g. spin 9 on chain I 
and 14 on II) and interplane A', for spins on chains in consecutive planes (e.g. spin 9 on chain I and 1 on III). 



Fig. 4. Magnetic excitations for the charge ordered (CO) phase, along 8 paths throughout the first Brillouin zone. F is taken 
as the energy unit. Parameters: F' = F — 1; A = 1 = A'; Si = 1.3, S2 = 1; c/ao — 1. 



Fig. 5. Magnetic excitations for the dimer phase, along 8 paths throughout the first Brillouin zone. Parameters: F 1 = 2, 
F = 1; A = 1 = A'; Si = S 2 = 1; c/a = L 
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Fig. 6. Magnetic excitations for the dimers phase, along 8 paths throughout the first Brillouin zone. For parameters: F' = 10, 
F = 1; A = 1 = A'; Si = S 2 = 1; c/a = 1. 
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